Neutron injection during primordial nucleosynthesis alleviates the primordial 7 Li problem 
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We present a parametrized study of the effects of free thermal neutron injection on primordial nucleosyn- 
thesis, where both the rate and the time scale of injection are varied. This generic approach is found to yield 
a successful solution for reducing the 7 Li abundance without causing significant problems to other elemental 
abundances. Our analysis demonstrates that hadronic injection, possibly due to decays or annihilations of dark 
matter particles with a mass of about 1 to 30 GeV, provides a possible solution to an outstanding problem in the 
standard Big Bang model. 



I. INTRODUCTION 

The motivation for this study concerns the discrepancy be- 
tween the primordial Li abundance predicted in the canon- 
ical Big Bang model and observational data. The primor- 
dial lithium abundance is deduced from observations of low 
metallicity stars in the halo of our Galaxy where the lithium 
abundance is almost independent of metallicity, displaying a 
plateau, the so-called Spite plateau |Q]]. This interpretation 
assumes that lithium has not been depleted at the surface of 
these stars, so that the presently observed abundance is sup- 
posed to be equal to the initial value. The small scatter of 
values around the Spite plateau is an indication that depletion 
may not have been very effective. Astronomical observations 
of these metal-poor halo stars H have led to a relative primor- 
dial abundance of: 

Li/H = (1.23+^) xlfj- 10 
A more recent analysis by Sbordone et al. H gives: 

Li/H = (1.58 ±0.31) x 1(T 10 . 
More generally, Spite and Spite 0] have reviewed recent Li 
observations and their different astrophysical aspects. Also 
see Frebel and Norris for a comprehensive review. 

On the other hand, the most recent Standard Big-Bang Nu- 
cleosynthesis (SBBN) calculations, using the most up-to-date 
nuclear data, give: 

Li/H = (5. 14 ±0.50) x 10 10 |H]. 
Hence there is a factor of 3-4 discrepancy between observation 
and theory at the WMAP7 baryonic density. 

7 Li is produced as a by-product of decay of 7 Be. Nuclear 
mechanisms to destroy this 7 Be have been explored. A possi- 
bly increased 7 Be(d,p)2oc cross section has been proposed by 
Coc et al. 10] and later by Cyburt and Pospelov (H] but was 
not confirmed by experiments JH-Qj]]. Other 7 Be destruction 
channels have recently been proposed by Chakraborty et al. 
lfl2tl and await experimental investigation. 

Another scenario would be to take advantage of an in- 
creased late-time neutron abundance, as introduced in ifllll for 
the generic case of hadronic injection. In the context of vary- 
ing constants, when the 'H(n,y) 2 H rate is decreased, the neu- 



tron late-time abundance is increased (with no effect on 4 He) 
so that more 7 Be is destroyed by 7 Be(n,p) 7 Li(p,oc)oc, [see in 
[Til Fig. 1]. Many other nuclear reactions could be poten- 
tial sources of free neutrons. However, a recent study lfl5ll 
extended the SBBN network to 59 nuclides from neutrons to 
23 Na, linked by 391 reactions involving n, p, d, t and 3 He in- 
duced reactions and 33 p-decay processes. The 7 Li abundance 
is now estimated to Li/H = 5.24 x 10~ 10 fiUl . as found also 
by llql . This confirms and even increases the discrepancy. 



Including physics beyond the standard model of particle 
physics and beyond the standard Big-Bang picture can also 
give rise to extra neutron injection. Indeed, BBN can be used 
as an anchor to test the plausibility of new physics, and con- 
versely, new physics can provide mechanisms to help solving 
the SBBN discrepancies with observations lfl7l - l20ll . One such 
option is that of hadronic decays of exotic unstable particles. 
For example, a metastable stop Next-to-Lightest Supersym- 
metric Particle (NLSP) decays into a gravitino Lightest Super- 
symmetric Particle (LSP), thus a dark matter candidate, and a 
top quark injects energetic protons and neutrons during nucle- 
osynthesis IIO, I2II - I27I1 . Another possible source of neutrons 
arises from residual annihilations of dark matter particles - 
such as neutralino LSP annihilating into fermion-antifermion 
couples that further hadronize - that are chemically decou- 
pled at BBN times 01 III III. In all these scenarios, neutron 
injection provides the primary impact on BBN and Li produc- 
tion. 



In this work, we study the effect of free neutron injection, 
parametrized by the injection rate and time-scale. Different 
injection models are thus included in the full code presented 
in Il5ll . Hence, this implementation is expected to give hints 
regarding the injection mechanism including possible nuclear 
reaction uncertainties, fundamental constant variations and 
exotic particle decays or annihilations. We comment on pos- 
sible scenarios behind neutron injection, however, we do not 
include a full treatment of the production and thermalization 
of neutrons in the code. 
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H. IMPLEMENTATION IN THE BBN CODE 

This code fl5ll is based on the Big Bang model and ACDM 
cosmology. There are three pieces of evidence for this physi- 
cal model: the universal expansion, cosmic microwave back- 
ground (CMB) radiation, and BBN. The latter comes from the 
primordial abundances of the "light cosmological elements": 
4 He, D, 3 He and 7 Li. They are produced during the first ~ 20 
minutes of the Universe when it was dense and hot enough for 
nuclear reactions to take place. The number of free param- 
eters entering the standard the BBN scenario has decreased 
with time. The number of light neutrino families is known 
from the measurement of the Z boson width by LEP exper- 
iments at CEPvN: N v = 2.9840 ± 0.0082 JH. The lifetime 
of the neutron enters in weak reaction rate calculations and 
many nuclear reaction rates have been measured in nuclear 
physics laboratories. The last parameter to have been inde- 
pendently determined is the baryonic density of the Universe. 
It is now deduced from the observations of the anisotropics 
of the CMB radiation coming from the Wilkinson Microwave 
Anisotropy Probe (WMAP) satellite. The number of baryons 
per photon, which remains constant during the expansion, 
r| is directly related to D.b by £V? 2 = 3.65 x 10 7 r|. The 
WMAP 7 year results now give Q. b h 2 = 0.02249 ± 0.00056 
andri = (6.16 ± 0.15) x 10 10 JH]. In this context, primordial 
nucleosynthesis is a parameter-free theory and is the earliest 
probe of the Universe. 

The main difficulty of the BBN calculations up to CNO is 
that an extensive network of reactions is needed, including «-, 
p-, a-, but also d-, ?-, and 3 He-induced reactions. Most of the 
corresponding cross sections cannot be extracted from exper- 
imental data only. In the BBN code, we use a more reliable 
rate estimates provided by the TALYS reaction code l30ll . 59 
nuclides are included, from neutron to 23 Na, linked by 391 re- 
actions involving «-, p-, d-, t-, and 3 He-induced reactions and 
33 (3-decay processes. 

Including an additional neutron injection in our SBBN code 
is straightforward. We allow protons to decay to neutrons with 
a lifetime of A. -1 (t ). As we are considering very low injection 
rates, this has no consequence on the high proton abundance. 
To illustrate the consequences of early or late injection, we 
consider the following cases: 

1. A(f ) = Ao at all time t (or temperature T) 

2. A(f ) = A for T < T c and for T > T c 

3. A(f ) = for T < T c and A for T > T c 

4. A,(f) = A, exp(-f/T v ) 

5. X(r) = X (£) 3 

with Ao constant and T c = 0.2 and 0.3 GK. Since the pro- 
ton abundance remains essentially constant (Y p «0.5 to 0.7) 
during BBN the rate of injection Y p (t)X(t) is constant in the 
(1), (2) and (3) intervals. Cases (4) and (5) represent more 
physical situations where neutrons come from the decay of a 
hypothetical particle X of lifetime T v decaying to X— >n + ... 
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FIG. 1: 4 He, D, 3 He and 7 Li abundances as a function of neutron 
injection rate for cases (1) (solid line), (2) with T c = 0.2 (sparse dots) 
and 0.3 GK (dots) and case (3) with T c = 0.2 (dash) and 0.3 GK 
(dash-dot). Hatched zones represent the observational outcome (see 
the text for details). 



with a branching ratio B„, or as a product of the annihilation 
of dark matter particles (discuss the relevant mass range of 
X is discussed below). In the latter cases the injection con- 
stant can be expressed as A-o = Y x (t = Q)B„/Y p z x . A cutoff 
in the neutron injection spectrum at lower temperatures (red- 
shifts) is expected due to the fact at a certain redshift (T ~ 
0.1-0.3 GK) the average time between interaction of neutrons 
becomes greater than the decay time of a neutron. 



III. RESULTS 

We summarize the outcome of the code in the figures. Fig- 
ureQ] shows that injection of neutrons at a rate of Ao ~ 10 -8 
s _1 when T > 0.3 GK alleviates the 7 Li problem without sig- 
nificantly affecting the other isotopes. Figure [2] shows that 
the 9 Be and n B abundances depend strongly on Ao while the 
CNO abundances are not modified. There is a modest (less 
than 50%) enhancement of D but this is well within astra- 
tion uncertainties. As one might expect, exotic particle decays 
with A = (1.5-2) x 10~ 8 s" 1 or (1-30) GeV dark matter an- 
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FIG. 2: Same as Fig.[I]but for 6 Li, 9 Be, 10 B, n B and CNO isotopes. 

nihilations with Xq = (3 — 5) x 10~ 9 s (Figure[3]l help solve 
the 7 Li problem. Figure @] shows the time evolution of the 
lightest elements for the SBBN and neutron injection from 
decaying exotic particle scenarios. The major impact of the 
injected neutrons is on 7 Be and 7 Li helping diminish the pri- 
mordial yield in 7 Li while the deviation on other light species 
show minor changes. 

In figures[T]and[3]the different observational constraints are 
given (green hatched zones). As previously stated, the Li/H 
abundance is obtained from 0]. The abundance of D is mea- 
sured in quasar absorption systems. The weighted mean value 
of the observations is D/H=(3.02 ± 0.23) x 1(T 5 (see JH 
for details). Note that two recent observations of D/H could 
slightly modify this value [32] and 13311 '). Finally, the deter- 
mination of the 4 He abundance in extragalactic H2 regions 
is fraught with difficulties due to systematic errors. Conse- 
quently, as shown in figures Q] and [3] the weighted mean value 
is Y p = 0.2566 ± 0.0028 still carries a large uncertainty 

IV. DISCUSSION 

While in our SBBN code, the neutrons are injected at 
equilibrium, it is likely that extra neutrons from any kind 
of beyond the Standard Model physics are produced out-of- 
equilibrium. It is, therefore, important to consider the ther- 
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FIG. 3: 4 He, D, 3 He and 7 Li abundances as a function of neutron 
injection rate for case (4) i.e. decay, with t v = 40 mn (solid) and 
case (5) i.e. annihilation with T c = 0.3 GK (dash). Hatched zones 
represent the observational outcome (see the text for details). 

malization process of neutrons during BBN before they decay, 
through which channels they do so, and to estimate the pos- 
sible perturbations to SBBN abundances. As mentioned by 
Jedamzik I23I l24ll . the thermalization process calculation of 
energetic nucleons is simplified greatly by two facts: first, the 
Hubble time is much greater than the mean time between any 

of the interactions under consideration Th rj 300 ( 9Q ^ eV ) " s, 
and second, the interactions between non-thermal and ther- 
mal nucleons are unlikely. The decay time of free neutrons 
(To = 881 s) is even greater than the Hubble time. 

There are 3 main classes of reactions: (1) elastic and inelas- 
tic n — p scattering, (2) the afore-mentioned spallation of 4 He 
with production of 3 He, and (3) both elastic and inelastic scat- 
tering n — 4 He. All of these processes contribute to thermal- 
ization, but the spallation to non-thermal 3 He might disturb 
the abundance of 6 Li 13511 through the following reactions 

n + 4 He -> 3 He + 2n 

4 He + 3 He -> 6 Li + p. W 

Firstly, we justify the claim that the injected neutrons in- 
deed thermalize before they decay and secondly, we estimate 
the production of 3 He and 6 Li. 
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f2 B h 2 =WMAP; X =10" 8 s" 1 ; T x =40 mn 




FIG. 4: 4 He, D, 3 He and 7 Li abundances in mass fraction as a func- 
tion of time for Xq = (dash) and case (4) with Xq = 10~ 8 s _1 
and X x = 40 mn (solid). Note that, in this case, 7 Li is lowered at 
X( 7 Li)= 6 x 1(T 10 which corresponds to Li/H= 1.1 x 10~ 10 . 



6 Li is a very interesting isotope as a new cosmological nu- 
cleus. Indeed, its abundance is measured in low metallicity 
stars and offers a unique probe of two different mechanisms of 
nucleosynthesis: SBBN and cosmic rays. The former produc- 
ing predominantly 7 Li, 6 Li was until recently considered as a 
pure spallative -i.e., post BBN- product. However, accord- 
ing to recent detections in very low metallicity PopII stars, 
the average is 6 Li/ 7 Li~ 0.042 13611 . These observations have 
been interpreted as evidence for a large primitive abundance 
of 6 Li ( 6 Li/H ~ 10~ 12 ) while SBBN calculations confirm a 
low primordial value ( 6 Li/H ~ 10~ 14 ). For details on the 
subject, see lfl5i l37ll . Recently, new studies using 3D atmo- 
sphere model in metal poor halo stars reconsider the detec- 
tion of 6 Li. In lT38h two detections are confirmed ( 6 Li/ 7 Li 
~ 5 — 10%). More observations are presently needed to im- 
prove the statistics. Nevertheless, this new spectroscopic re- 
search can be an indicator of new physics, as has been point 
out by Jedamzik l39ll . In that regard, we have to make sure 
that the extra physics we consider does not perturb signifi- 
cantly the abundance of 6 Li. We do it in an order of magnitude 
estimate. Figure |2] shows clearly that our models should not 
modify too much 6 Li: at Ao — 10~ 8 n s _1 6 Li is in the range 



The total cross-section of elastic and inelastic scatterings of 
n off p is about 70 mb at 100 MeV and 30-40 mb above 100 
MeV. The cross-section of the reaction n + 4 He — > 3 He + 2n 
varies from about 15 mb at 30 MeV and 50 mb at 50 MeV 
down to 20 mb at 1 GeV. The sum of the cross-sections of 
the remaining inelastic processes is comparable to n + 4 He — > 
3 He + 2n at 50 MeV and about 3-4 times greater than the rates 
above 100 MeV. The elastic scattering of neutrons off 4 He 
varies from 500 mb at 20 MeV down to 30 mb at 100 MeV 
and stays constant at higher energies. 

The characteristic time of scattering of a neutron of kinetic 
energy Ek off thermal H and 4 He is 



i A, 



\ -3 



( 90 KeV J ^50mb 



-l 



l + (m n /50MeV) 2 
1 + (m n /E k ) 2 



-1/2 



s, 



where Ah = 26 and A* He = 2.17, is the kinetic energy 
of the neutron. In the relativistic limit, E„ ^> m n , x goes down 
to 0.07 s, while the decay time grows linearly with energy 

We also have to make sure that a significant fraction of en- 
ergy is transferred in a single scattering or per mean free path 
length (its inverse is denoted E (XdE / dx) ~ 1 in Jedamzik 12411 ). 
In inelastic n — p scatterings, ^ remains constant up to 250 
MeV (^f > 0.1). Indeed in elastic n - 4 He scatterings, the 
average energy transfer is about 5 — 10 MeV. In inelastic pro- 
cesses n — 4 He, the neutron loses as much as the binding en- 
ergy of 4 He (28.3 MeV) and a quarter of the remaining energy. 

From these arguments, it is clear that extra neutrons ther- 
malize in these conditions before decaying in a range of tem- 
peratures from 100 KeV down to a few KeV and in the range 
of energies from about 10 MeV to 1 GeV. 

The hypothesised extra neutrons might be produced by an- 
nihilating dark matter particles. The energy injection due to 
dark matter annihilations if the freeze-out of the dark matter 
species happened at BBN temperature is severely constrained 
(see lfl8ll for example). However, if the freeze-out happened 
before BBN, annihilations become marginal, as the expansion 
rate dominates the interaction rate. Nonetheless there would 
be a residual annihilation rate of dark matter into standard 
model particles. Eventually, after hadronization, a spectrum 
of neutrons would be generated, that would reach thermal 
equilibrium as discussed before. 

The annihilation rate of uniformly distributed dark matter 
per baryon can be written as 



Tb = i(ov)^M(l+z) J 



5.3 x 10- 



(ov> 



3xl(T 



30 GeV V 
M X ) 



^90keV 



)V 



(1-6) x 10 



-14 



(2) 

where iiq.x and «o j, are the present day number densities of 
dark matter and baryons. We see that at a temperature of about 
90 keV, a particle dark matter mass of Mx r$ 30 GeV and a 
canonical annihilation rate are plausible parameters needed 
to achieve X = (3 — 5) x 10~ 9 s _1 , depending on the neutron 
spectrum generated by the annihilations. 

The neutron spectrum is generated after hadronization of 
the particles produced at annihilation, and it is expected to be 
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peaked at roughly Mx/5 — Mx/15. Therefore, the lighter the 
dark matter particle, the larger the fraction of thermal neu- 
trons. However, the dark matter has to be heavy enough to 
produce neutrons, hence, the most interesting mass range lies 
roughly between 1 and 30 GeV. 

The annihilation rate and branching ratios depend on the 
dark matter candidate. Moreover, the dark matter tempera- 
ture evolves from chemical decoupling down to thermal de- 
coupling (see l40|p . The dependence of the annihilation rate 
on the dark matter temperature can be very strong; for ex- 
ample, if the freeze-out mechanism invokes a nearly resonant 
exchange, or co-annihilations 14111 . 

A relevant example for a dark matter candidate in the mass 
range discussed here is the neutralino in the Next-to-Minimal 
Supersymmetric Standard Model. As shown in ll42ll . the res- 
onant mechanism at freeze-out can yield a very large boost to 
the annihilation rate at lower temperatures (see Fig. 4 in |42i]). 
For kinetic decoupling at ~ 7y / '10, one could have a fac- 
tor 10—100 enhancement in the annihilation rate from the 
3 x 10 _26 cm 3 s _1 required at freeze-out. Thus there can be a 
variety of dark matter candidates (with different masses and 
annihilation cross-section mechanisms) which provide an in- 
jection flux of (3 —5) x 10~ 9 s -1 . 

It is interesting to note that some of these candidates could 



explain direct detection signals as they have the right mass 
range and could attain the needed interaction rates with de- 
tectors. Also, they could be challenged by y-ray production 
at dwarf spheroidal galaxies. Relating annihilating rates at 
freeze-out, BBN and galactic times, and elastic scattering in- 
teractions with nuclei, can provide powerful constraints on a 
given dark matter model. 

In conclusion, neutron injection can help to resolve the 7 Li 
problem provided that the neutrons are essentially thermal. 
This can be achieved for annihilations or decays of dark mat- 
ter particles in the mass range 1-30 GeV. A detailed physi- 
cal model involving, for example, a metastable supersymmet- 
ric NLSP or annihilating neutralino dark matter is beyond the 
scope of this paper, but would seem to be easily achievable. 
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